Length-dependent conductance and thermopower in single-molecule junctions of 

dithiolated oligophenylene derivatives 
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We study theoretically the length dependence of both conductance and thermopower in metal- 
molecule-metal junctions made up of dithiolated oligophenylenes contacted to gold electrodes. We 
find that while the conductance decays exponentially with increasing molecular length, the ther- 
mopower increases linearly as suggested by recent experiments. We also analyze how these transport 
properties can be tuned with methyl side groups. Our results can be explained by considering the 
level shifts due to their electron-donating character as well as the tilt-angle dependence of conduc- 
tance and thermopower. Qualitative features of the substituent effects in our density-functional 
calculations are explained using a tight-binding model. In addition, we observe symmetry-related 
even-odd transmission channel degeneracies as a function of molecular length. 

PACS numbers: 85.65.-|-h, 65.80.-l-n, 73.23.Ad, 73.63.Rt 



I. INTRODUCTION 

In the field of molecular electronics, research has so far 
mostly concentrated on the dc electrical conduction prop- 
erties of single- molecule contacts.^ By now it is known 
that the charge transport through organic molecules is 
typically due to electron tunneling. This is evidenced in 
particular by the exponential decay of the conductance G 
with increasing molecular length in contacts formed from 
oligomers with varying numbers of unitsi^i^ii^ Consider- 
ing the statistical nature of experiments at the molecu- 
lar scale, a conclusive comparison to theory is presently 
difficult. Nevertheless, associated decay coefficients ap- 
pear to be reproduced by theoretical calculations >^i^i2, For 
a deeper understanding of molecular-size contacts, it is 
useful to analyze also other observables in parallel with 
the dc conductance. Emerging new lines of research in- 
volve the photoconductance^iSiiiS and heat transport.— 
In this paper we concentrate on the thermopower. For 
metallic atomic contacts, this quantity was studied ex- 
perimentally already some years ago,-'' but for molecular 
contacts only very recentlyJ^ii^ 

The thermopower Q, also known as the Seebeck co- 
efficient, measures the voltage AV induced over a con- 
ducting material at vanishing steady state electric cur- 
rent /, when a small temperature difference AT is ap- 
plied: Q = -AV/AT\j^Q. It is known that in bulk 
materials the sign of the thermopower is a hint of the 
sign of the main charge carriers. If Q < {Q > 0), 
charge is carried by electron-like (hole-like) quasiparticle 
excitations, as in an n-doped (p-doped) semiconductor.^® 
Analogously, the thermopower of the molecular junction 
gives information about the alignment of the energies 
of the highest occupied and lowest unoccupied molec- 
ular orbitals (HOMO and LUMO) with respect to the 
metal's Fermi energy EpJ^-^^^ In the experiment,^** Q 



was measured for gold electrodes bridged by dithiolated 
oligophenylene molecules. It was found to be positive, 
which indicates that Ep lies closer to the HOMO than to 
the LUMO. Also, it was observed that Q grows roughly 
linearly with the number TV of the phenyl rings in the 
molecule. More recently*^ the effects of substituents and 
varied end groups have been analyzed for the benzene 
molecule. 

Using transport calculations based on density- 
functional theory (DFT) we investigate in this paper 
the length dependence of the conductance and the See- 
beck coefficient for dithiolated oligophenylenes bonded 
to gold contacts. For the molecules studied in Ref. IT^ . 
we find that the conductance decays exponentially with 
increasing molecular length. Decay coefficients compare 
reasonably with previous DFT calculations""'®'^ and also 
with experiments)^!^ considering differences in contact 
configurations, molecular end groups, and uncertainties 
with respect to environmental effects. In addition, the 
thermopower increases linearly, with a magnitude co- 
inciding with the measurements We also study how 
the results change as various numbers of methyl sub- 
stituents are introduced to the molecules. The effect 
of these substituents is twofold: (i) they push the en- 
ergies of the IT electrons up as a result of their electron- 
donating behavio r^^'^" and (ii) they increase the tilt an- 
gles between the phenyl rings through steric repulsion. 
The latter effect tends to decrease both G and Q due 
to a reduction of the degree of 7r-electron delocalization, 
while the former opposes this tendency by bringing the 
HOMO closer to Ep- A simplified 7r-orbital model is used 
to explain essential features of the DFT results. 

The paper is organized as follows. Sec. [Til explains 
details of our DFT-based approach. Then, in Sec. lIIIl we 
present the molecular contacts, whose charge transport 
properties we determine in Sec. IIVI In Sec. |V] we show. 
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Figure 1: (Color online) The molecules studied in this work. 
When they are contacted to the gold electrodes, sulfur atoms 
replace the terminal hydrogen atoms. 

how the DFT results can be understood in terms of the 
TT-orbital model, but we discuss also effects beyond this 
simplified picture in Sec. IVII Finally, we end in Sec. IVIII 
with a discussion and conclusions. 



II. METHODS 

The general formulas for the treatment of thermoelec- 
tric effects, based on the Landauer-Biittiker formalism, 
are discussed in detail in several references] ^^i^^'^^i^^'^^i^^ 
By expanding the expression for the current / to linear 
order in AV and AT and considering the cases AT = 
and / = 0, respectively, one arrives at 

G^GoKoiT), Q = -^}fl- (1) 

with K^iT) = jT{E)[~df{E,T)/dE]dE, K,{T) = 
J {E - fi)T{E)[^df{E,T)/dE]dE, and Gq = 2e^/h. 
Here t{E) is the transmission function, f{E,T) = 

{exp [(i? — /i) /fc^T] + 1} ^ the Fermi function, and /i 
the chemical potential, w Ep. At low temperature, the 
leading-order terms in the Sommerfcld expansions yield 

G^GoriE,), Q = -^^^fcBT, (2) 

e 3 t[Ef) 

where prime denotes a derivative. In our DFT-based re- 
sults presented below, we calculate G and Q according 
to Eq. ([1]). However, Eq. ([2]) approximates the results 
to within a few percent at room temperature, since t{E) 
is smooth around Ep due to the off-resonant situation. 
It should be noted that the equations neglect electron- 
vibration interactions. These could in principle be in- 
cluded by adding the inelastic corrections to the expres- 
sion for the current jSi' but we expect also these contribu- 
tions to be relatively small even at room temperature. 

The transmission functions are computed with the help 
of Green's function techniques. The electronic structure 
is described in terms of DFT as implemented in the quan- 
tum chemistry program TuRBOMOLE V5.7, where we 
employ the BP86 exchange-correlation functional and the 
standard Gaussian basis set of split valence quality with 
polarization functions on all non-hydro gen atomsi^ For 
further details on our method, see Refs . I26ll27ll28ll29l . 
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Figure 2: (Color online) HOMO and LUMO energies of the 
molecules of Fig. [1] 



III. CONTACTS 

The molecules studied are shown in Fig. [T] Those 
labeled with Rl to R4 are the pure oligophenylenes, 
S2 to S4 denote oligophenylenes where the hydrogen 
atom in one of the two ortho positions with respect to 
each ring-connecting carbon atom is substituted with a 
methyl group^Si^i and D2 to D4 have substituents in 
both ortho positions. Here, the numbers N = 1,...,4 
refer to the number of phenyl rings in the molecule. 
The tilt angles for the R, S, and D molecules vary be- 
tween 33.4° < ipR < 36.4°, 84.8° < (fs < 90°, and 
88.8° < fD < 90°, and the distances between the ter- 
minal carbon atoms of the molecules are described to a 
good accuracy by d = a -I- bN, with a — —0.154 nm and 
b ^ 0.435 nm. 

The HOMO and LUMO energies of the molecules are 
shown in Fig. ^ It can be noticed that the HOMO- 
LUMO gaps of the S and D series are larger than those 
of the R series. 

To form the junctions, each molecule is coupled to the 
hollow position of the tips of two gold [111] pyramids via 
a sulfur atom. The atomic positions of the molecule and 
the first gold layers of the tips are then relaxed. This is 
depicted in Fig. [3] for S3. The relaxed part is also the 
"central region" in the transport calculationsi^S The tilt 
angles ip and the distances d of the contacted molecules 
are not essentially different from those of the isolated 
molecules. 




Figure 3: (Color online) (a) The contacts to the gold elec- 
trodes are formed through sulfur atoms bonded to the hollow 
position of the tip of a [111] pyramid, (b) For the transport 
calculations, the tips of the pyramids are part of the extended 
molecule (C). The more remote parts (blue) are absorbed into 
ideal semi-infinite left {L) and (R) surfaces (gray). 
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Figure 5: (Color online) The parameters of the TB model 
represented by a molecule with A'^ = 3 phenyl rings. Also 
shown are schematic graphs of the onsite energies ej of the 
different rings j — 1, . . . ,N for R, S, and D type molecules 
with A'' = 1, . . . , 4 rings. 



Figure 4; (Color online) (a,b) Transmission function and the 
negative of its logarithmic derivative. (c,d) The corresponding 
G and Q, including the temperature corrections in Eq. ([ij. 
The experimental data in (d) are from Ref. Ilil . The straight 
lines are best fits to the numerical results for R (solid line), S 
(dashed line), and D (dash-dot-dotted line). 



IV. DENSITY-FUNCTIONAL-BASED 
TRANSPORT 

The results for the transmission, its logarithmic deriva- 
tive, the conductance, and the thcrmopower for all of 
the R, S, and D type molecular junctions are collected in 
Fig. m In order to compare with the room-temperature 
experiments^^ we set T = 298 K in Eq. ^ . The Seebeck 
coefficients are displayed in Fig. IH^d) together with the 
experimental results, where the molecules Rl, R2, and 
R3 were studied. 

With increasing number N of phenyl rings, the con- 
ductance decays as G/Go ^ e^^^ [Fig. [DJc)]. For the 
R series, we find the decay coefficient Pr = 1.22f^ in 
good agreement with theory of Ref. 0.— Previous the- 
oretical estimates^ and experimental results^i^ for thi- 
olated self-assembled monolayers are consistently some- 
what larger than this value. In particular, the results 
reported in Ref. [3 vary between 1.5 < Pr < 2.1, and 
for amine end groups an experimental value of fin = 1.5 
was reported)^ The underestimation of may indicate 
an overestimation of the conductance computed within 
DFT^i24 However, a comparison between theory and 
experiment is complicated due to the differences in the 
end groups usedi^ Therefore, no conclusive comparison 
is possible. For the S and D series, we find Ps = 3.69 
and Pd — 4.07, which are both much larger than Pa. 
This increase reflects the reduced delocalization of the 
TT-electron system. The absolute conductance values of S 
and D are very similar to each other. 

Since the Fermi energy at Ep = —5.0 eV lies closer to 
the HOMO than to the LUMO level, the Seebeck coeffi- 
cient Q has a positive value [Fig.lD^d)]. We also find that 
Q increases roughly linearly with N, as suggested by the 
experiments.^'* Indeed, assuming that the transmission 



around E = Ep is of the form t{E] = a{E)e-'^'^^^^ , 
then Eq. © yields Q = Q(°) + Q'-^'N, where Q^^^ = 
-klTTT^[\naiEF)]'/3e and Q(i) = klTn^ P' {Ep) /ie.^^ 
Two things should be noted here. First, Q does not nec- 
essarily extrapolate to zero for TV = 0, leading to a fi- 
nite "contact thermopower" Second, depends 
on the prefactor a{E), but Q^^^ does not. Since a{E) 
contains the most significant uncertainties related to the 
contact geometries, Q^^^ can be expected to be described 
at a higher level of confidence than Q'"'- Best fits to 
our results and the experimental data give Q^'' = —0.28, 

Q^^^ = 7.77 and Q^.^ = 6.43, Q^.^ = 2.75 ^V/K, 
respectively. Differences in the fit parameters mainly 
stem from the data point for R3, where the experimen- 
tal value is lower than calculated. Considering the re- 
ported order-of-magnitude discrepancies between mea- 
sured conductances and those computed from DFTj^i^i 
the agreement still appears reasonablei^ However, for 
large enough T and the abovementioned exponential 
and linear laws for the length dependences of G and Q 
should be modified due to interactions with the thermal 
environment The fact that the experimental data in 
Fig. mjd) exhibit a rather good linearity suggests that 
molecular vibrations do not play a crucial role and that 
the electronic contribution to Q is dominant. 

Although the conductances for S and D are very sim- 
ilar, their thermopowers are rather different. Further- 
more, the magnitudes for N > 2 follow the surprising 
order Qs < Qr < Qd-^^ As we will discuss below, these 
observations can be understood through the two com- 
peting substituent effects: (i) a change in the alignment 
of the TT-electron levels with respect to Ep as a result 
of the electron-donating nature of the methyl group ^^i^° 
and (ii) an increase in the ring-tilt angles. For the iso- 
lated molecules the second effect results in the opening 
of the HOMO-LUMO gap when going from molecules R 
to S or D, while the first effect causes the difference be- 
tween the HOMO and LUMO energies of the S and D 
series (Fig. [5]). 
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molecule HOMO-1 



:^ i|t 3^ 

-6.167 -6.265 eV -1.125 eV -1.124eV 

J*Y^ ^^pk J^4ft. M^s. 

-6.174 eV -5.948 eV -1.093 eV -0.982 eV 



benzene 



toluene 




meta-xylene -5.969 eV -5.779 eV -1.016 eV -0.906 eV 







para-xylene -6.088 eV -5.678 eV -0.987 eV -0.938 eV 




durene -5.756 eV -5.393 eV -0.842 eV -0.699 eV 

Figure 6: (Color online) Influence of methyl substituents on 
the level alignment of frontier molecular orbitals in benzene. 
From left to right the molecular structure is displayed together 
with isosurface plots of the HOMO-1, HOMO, LUMO, and 
LUMO+l wave functions. Below these plots the common 
name of the molecules and the energies of the molecular levels 
are indicated. 



V. TT-ORBITAL MODEL 

In order to gain some understanding of the general 
features of the dependence of G and Q on the num- 
ber of substituents, we study a simple tight-binding 
(TB) model, which describes the 7r-electron system of 
the oligophenylenes (Fig. [5]). The onsite energies ej 
(j = 1, . . . , N) are equal on all carbon atoms of phenyl 
ring J, the intra- ring hopping t is assumed to be the same 
everywhere, and the inter-ring hopping u is parametrized 
through u — t cos If. We assume the effect of the side 
groups to come into play only through ip and cj. The 
leads are modeled by "wide-band" self-energies El^^ = 
— zF, acting on the terminal carbon atoms. 

We extract the parameters of our model as follows. 
For the R molecules, we set ej — 0. In the S and D 
molecules, the e^'s for rings with one, two, or four methyl 
groups are obtained by performing DFT calculations of 
methylbenzene, dimethylbenzene, and tetramethylben- 
zene (Fig. (H). We find that the HOMO and LUMO 
energies in these molecules shift upwards monotonously 
with the number of methyl substituents. This can be at- 
tributed to the electron-donating character of the methyl 
groups, and the resulting increase of Coulomb repulsion 
on the phenyl ring. We use the averages of the HOMO 
and LUMO shifts relative to benzene, and obtain, respec- 
tively, e^""^) = 0.17, e(*™°) = 0.33, and e(^°"'~) = 0.58 
eVj^ The onsite energies of the different phenyl rings 
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Figure 7: (Color online) (a) Transmission functions for the 
TB model of Fig. [5] with parameters chosen as explained in 
the text, (b) Corresponding Seebeck coefficients according to 
Eq. ([2]). The straight lines are best fits as in Fig. |4l (c,d) G 
and Q as a function of the tilt angle for the molecules with 
A'^ = 3. The vertical dotted lines indicate the "equilibrium 
angles" ipR, ips, and ipn, and the arrows represent changes 
when going from R3 to S3 to D3. 



of the R, S, and D molecules with N = 1,...,4 are 
schematically represented in the lower part of Fig. [S] The 
hopping t is set to half of (the negative of) the HOMO- 
LUMO gap of benzene, with the result t = —2.57 eV 
(Fig. [2|). The scattering rate F and the tilt angles for 
R, S, and D molecules are chosen to reproduce approx- 
imately the minimal transmissions for Rl, R2, S2, and 
D2 in the DFT results of Fig. H^a). This yields F = 0.64 
eV, ipB. = 40.0°, ips ^ 85.5°, and tpo = 86.5°. The last 
free parameter is Ep, which should be determined by the 
overall charge-transfer effects between the molecule and 
the electrodes. We choose its value of Ep = —1.08 eV 
close to the crossing points of the transmission curves of 
the S and D molecules [cf. Figs. Ufa) and[7I^a)]. 

Similar features can be recognized between the DFT 
results and our model. In particular, as shown in 
Fig. [Tjb), the correct order of the thermopowers is re- 
produced. In going from R to S, the increase in the 
tilt angle and the associated breaking of the 7r-electron 
conjugation dominates over other side group effects. As 
a result both G [Fig. [DJc)] and Q [Fig. Uld)] decrease. 
In going from S to D, an interplay of the sidegroup- 
induced level shifts and a small residual increase in tp 
raises Q above the value for R, while G remains al- 
most unchanged. We note that for an A^-ring junction 
the lowest-order terms in an expansion of the transmis- 
sion in powers of u yield (/^-dependences of the form 
t{E) « ci{E) cos2(^-i) ip + C2{E) cos^^ if, for E « Ep.^^ 
Independently of N, this results in Q qi + q2 cos^ (p, 
where we find qi,q2 > in our case. For N = 2, the 
neglect of the C2-term leads to the well-known "cos^ cp 
law" of Gr^ but we see that for the tilt-angle depen- 
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Figure 8: (Color online) (a) Tilt-angle-dependent transmis- 
sion T = Ti resolved in its transmission channels Ti for the 
molecule R2 and (b) thermopower for R2, S2, D2. The curve 
for R2 is fitted by a function of the form gi + q2 cos^ cp with 
gi — 13.27 and q2 — 1.38 /iV/K for the tilt-angle interval 
from 0° to 60°. Dotted vertical lines indicate DFT equilib- 
rium tilt-angles ipR2, 'fiS2, and ipD2- 



dencG of Q the presence of this term is significant, since 
(72 = if C2 = 0. When Ep is close to a resonance, fur- 
ther higher-order terms become increasingly important 
and deviations from the qi + q2 cos^ ip law result. This is 
seen most clearly as the non-monotonous iy9-dependence 
of Q for molecule D3 in Fig. [ltd) . 



VI. EFFECTS BEYOND THE tt-ORBITAL 
MODEL 

Close to perpendicular ring tilts results from the tt- 
orbital model should be taken with care. At Lp = 90° 
the TT-TT coupling u vanishes, and any other couplings 
between the rings will become important.'**' Let us an- 
alyze this for the biphenyl molecules R2, S2, and D2. 
For them, we have varied the tilt angle between the 
rings, and have determined the charge transport prop- 
erties for every Lp using our DFT-based approach. Fur- 
ther details of the procedure are described in Ref. [2^. 
For R2 (and similarly for S2 and D2) we observe that 
for most tilt angles [ip < 80°) the transmission is dom- 
inated by a single channel [Fig. [Hl^a)] of tt-tt character. 
However, at large angles two transmission channels of 
the same magnitude are observed, which become degen- 
erate at 1^ = 90°. They arise from ct-tt couplings between 
the two rings, whose strengths are proportional to sini^. 
Thus the two degenerate channels are of the ct-tt and tt-ct 
type. These features are obviously not accounted for by 
the TB model, where only a single tt-tt channel is present 
independently of ip. The ct-tt couplings should also mod- 
ify the tilt-angle dependence of the thermopower, which 
is plotted in Fig. HUb) . Due to the electron-donating na- 
ture of the methyl groups the thermopower at fixed Lp 
increases from R2 to S2 and D2. While the curve for R2 
can be described by the law gi -I- cos^ </3 for (/? < 80°, a 
dip is observed for larger (p. Similar deviations in Q are 
also present for S2 and D2, where we have investigated a 
smaller tilt-angle interval because of the steric repulsion 
of the methyl groupsi^ 
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Figure 9: (Color online) Importance t\/t2 of the second trans- 
mission channel T2 as compared to the first ri as a function 
of molecular length for the molecules of series R, S, and D. 



The degeneracy of the transmission channels is due to 
the D2d symmetry of biphenyl when tp — 90° 42. For the 
longer oligophenylenes the symmetry D2d can occur if 
and only if tilt angles are all at 90° and the number of 
rings is even. Hence for S2, S4, D2, and D4 the ratio 
T2/T1 of the first two transmission channels should be 
particularly large. Such even-odd oscillations are indeed 
visible in Fig. [5] for the S and D series, while they are 
absent for R. Owing to the fact that tilt-angles deviate 
from 90°, the oscillations decay. In particular T2/T1 is 
much smaller for S4 than for D4, where the minimal tilt 
angles in the contacts are 85.1° and 89.1°, respectively. 



VII. DISCUSSION AND CONCLUSIONS 

As is well known, there are many theoretical uncer- 
tainties involved in the determination of transport prop- 
erties based on DFT calculations, and improvements for 
the methods of molecular-scale transport theory are cur- 
rently being soug htS^iMAL Indeed, it is typical to find 
order-of-magnitude differences between measured con- 
ductances and those computed within DFTi^iM It is also 
known that atomic configurations can have a strong in- 
fluence on the conductance,— i^i^ and hence a conclu- 
sive comparison with experimental data would require 
the statistical analysis of a large number of contact ge- 
ometries. However, models predict Q to be insensitive 
to changes in the lead couplingSfii and hence it is ex- 
pected to be a more robust quantity that G. Although 
the measurements of Ref. |lj were carried out at room 
temperature, we do not expect molecular vibrations to 
play an essential role in the results due to the weakness 
of the electron- vibration coupling. Thus, the comparison 
we have made with our elastic transport theory seems 
justified. Nevertheless, the investigation of the effect of 
molecular vibrations is an interesting direction for fu- 
ture research, also in view of optimizing the properties 
of molecular thermoelectric devices for potential appli- 
cations. As far as gaining a basic understanding is con- 
cerned, low-temperature measurements would be desir- 
able to remove uncertainties related to the thermal exci- 
tation of vibrations. The purpose of the present work was 
to study the general trends for a series of molecules that 
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are similarly coupled to the electrodes, and considering 
all the potential uncertainties, the agreement we obtain 
for the thermopower appears to be quite reasonable. 

In conclusion, we have analyzed the length depen- 
dence of conductance and thermopower for oligopheny- 
lene single-molecule contacts. While we found the con- 
ductance to decay exponentially with the length of the 
molecule, we observe that the thermopower increases lin- 
early. For possible future applications it is interesting 
to know, how the magnitudes of these quantities can be 
tuned. We analyzed how this can be achieved by chemical 
substituents, which control the Fermi-level alignment and 
the molecular conformation. We demonstrated that a 
simple TT-electron tight-binding model can help to under- 
stand basic substituent effects. In addition, we observed 
an even-odd effect for transmission channel degeneracies 



upon variation of the number of phenyl rings, which we 
explained by molecular symmetries. 
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